Entanglement Across a Transition to Quantum Chaos 



in 
o 
o 

(N 



Carlos Meji'a-Monasterio^Q Guliano Benenti( a ' b < c )|j] Gabriel G. Carlo( a ' fc )@ and Giulio Casati( a - fc ' c J! 
Center for Nonlinear and Complex Systems, Universita degli Studi dell'Insubria, via Vallegio 11, Como 22100, 

Istituto Nazionale per la Fisica della Materia, 
Unitd di Como, via Vallegio 11, Como 22100, Italy and 
^Istituto Nazionale di Fisica Nucleare, Sezione di Milano, Via Celoria 16, 20133 Milano, Italy 

(Dated: February 1, 2008) 

We study the relation between entanglement and quantum chaos in one- and two-dimensional 
spin-1/2 lattice models, which exhibit mixing of the noninteracting eigenfunctions and transition 
from integrability to quantum chaos. Contrary to what occurs in a quantum phase transition, the 
onset of quantum chaos is not a property of the ground state but take place for any typical many-spin 
quantum state. We study bipartite and pairwise entanglement measures, namely the reduced Von 
Neumann entropy and the concurrence, and discuss quantum entanglement sharing. Our results 
suggest that the behavior of the entanglement is related to the mixing of the eigenfunctions rather 
than to the transition to chaos. 
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I. INTRODUCTION 

Quantum entanglement has been identified as a key 
ingredient in quantum communication and information 
processing. The content of entanglement in a particu- 
lar system is considered as a resource to perform several 
tasks in a more efficient and more secure way than any 
other classical method For instance, quantum tele- 

portation protocols require to share a maximally entan- 
gled state between the sender and the receiver. In the 
field of quantum computation, it has been found that for 
the case of quantum algorithms operating on pure states, 
the presence of multipartite entanglement between the 
components constituting a quantum processor is a nec- 
essary condition to achieve an exponential speedup over 
classical computation 

On the other hand, for the operability and stability of 
any quantum computer, the entanglement can also play 
the role of the property to be minimized. The unavoid- 
able entanglement between the quantum processor and 
the environment is one of the most important sources 
of noise and, therefore, of computational errors. The 
understanding and control of noise in quantum proto- 
cols is clearly needed to implement any reliable quantum 
computation 0, 0, IE IE @- Even when a quantum 
processor is ideally isolated from the environment, i.e., 
in situations where the decoherence time of the proces- 
sor is very large as compared to the computational time 
scales, the o per ability of the quantum computer is not yet 
guaranteed [l(| ■ Indeed, also the presence of device im- 
perfections hinders the implementation of any quantum 
computation. A quantum computer is a quantum many- 
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body system. The interaction between the qubits com- 
posing the quantum registers of the computer is needed 
to produce the necessary amount of entanglement. More- 
over, device imperfections like small inaccuracy in the 
coupling constants induce errors. Above a certain imper- 
fection strength threshold (chaos border), quantum chaos 
sets in [Tol [ill IT^L fl3| . In such a regime, exponentially 
many states of the computational basis are mixed after 
a chaotic time scale. This sets an upper time limit to 
the stability of a generic superposition of states coded 
in the quantum computer wave function. A necessary 
requirement for quantum computer operability and fault 
tolerant computation schemes is the possibility to oper- 
ate many quantum gates inside the chaotic time scale. 

For the case of many-body systems, the transition to 
chaos has been studied for fermions and bosons (see e.g. 
|l4| and references therein), and particularly for lattice 
spin systems [l5j |. In fermionic systems with two-body 
interactions it has been found that, if the interaction 
strength exceeds some critical value, fast transition to 
chaos occurs in the Hilbert space of many-particle states 
|16| . This is commonly studied in terms of the transition 
between the different spectral statistics of integrable and 
chaotic systems. For systems with a finite size, this tran- 
sition is smooth and only a crossover border where the 
transition occurs can be identified. The question whether 
this smooth transition becomes sharp in the thermody- 
namic limit is still under debate. However, in some cases 
a sharp transition to chaos is found, e.g., in the three- 
dimensional Anderson model [l7|. 

In a different context, the behavior of quantum en- 
tanglement across a quantum phase transition has re- 
cently attracted much attention. Quantum phase transi- 
tions (QPT) consist in a qualitative change in the ground 
state of the system as some control parameter is varied 
|l8j. Unlike classical phase transitions, QPT occur at 
zero temperature and the fluctuations developed at the 
critical point are fully quantum. These fluctuations dom- 
inate the behavior of the system near the critical point, 
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where correlations are long-range in character. It has 
been recently pointed out that the genuine quantum char- 
acter of quantum phase transitions is due to entangle- 
ment 0] . It has also been argued that the ground state 
of the system is strongly entangled at the critical point 
[Ifi, l2fT|. Therefore, the behavior of entanglement across 
QPT is particularly interesting for quantum computation 
and communication, where a maximization of the content 
of entanglement is desirable. The study of the relation 
between QPT and entanglement has been focused on the 
possible universal behavior of the entan glem ent content 
at the transition |H El El El El ElM El Hj. In 
this context, it has been shown that in different model 
systems quantities like the derivative of the concurrence 
[2TL 1221 ] and Von Neumann entropy [24l l3fT | present criti- 
cal behavior at the transition. The dependence of entan- 
glement on disorder and its interplay with chaos has also 
been studied [3ll I32I l33| . Moreover, the evolution of en- 
tanglement in quantum algorithms simulating quantum 
chaos has been recently investigated [3>E1IH1I3 

The aim of this paper is to characterize the behav- 
ior of quantum entanglement in non-integrable systems 
when a transition to quantum chaos occurs. We are in- 
terested in understanding how the entanglement content 
behaves in transitions from integrability to chaos. We 
would like to stress that, differently from the previous 
studies on QPT HH H3, 123, H3, HE EE E3, EM IS^I , the 
transition to chaos is not only a property of the ground 
state but takes place for any typical many-body state. 
We numerically study two lattice models of interacting 
many spins that show a transition to chaos. They have 
previously been studied as models of the quantum com- 
puter hardware [H El 111 

In both models, the 
transition to chaos is driven by the strength of the in- 
teraction between the spins. We consider bipartite and 
pairwise entanglement measures and focus on the rela- 
tion between the entanglement and the onset of chaos. 
We focus our study on the eigenstates corresponding to 
the center of the spectrum, where the many-body den- 
sity of states is larger and therefore quantum chaos sets 
in at small interaction strengths. Nevertheless, also the 
behavior of the other parts of the spectrum is discussed. 
We use exact diagonalization techniques to obtain all the 
eigenstates of the considered spin models. Therefore, we 
are limited to consider relatively small system sizes, from 
which the study of any possible finite-size scaling for the 
behavior of the entanglement at the chaos border is out 
of reach. Nevertheless, we discuss, at a qualitative level, 
the similarities and differences between the behavior of 
entanglement across a QPT and at the onset of quantum 
chaos. We show that the dependence of pairwise entan- 
glement on the size, distance and range of the interactions 
can be understood in terms of the sharing of entangle- 
ment among the different parties of the system. Further- 
more, we demonstrate that the behavior of entanglement 
is related to the mixing of noninteracting eigenfunctions 
rather than to the transition to chaos. 

This paper is organized as follows. In Sec. |H] we re- 



view the measures of entanglement that we will use. The 
measures that signal the onset of quantum chaos are re- 
viewed in Sec. IIIII In Sec. IIVI we define the spin lattice 
models investigated in this paper and discuss numerical 
results on their entanglement content and their transition 
to chaos. In Sec. [3 we present our final remarks. 

II. MEASURES OF ENTANGLEMENT 

A pure state is said to be separable if for a given 
partition of its Hilbert space H — Ha <8> Hb it can be 
written as \ip) = |o) ® \b). Here |o) and \b) are vectors re- 
siding in the Hilbert subspaces Ha and Hb, respectively. 
The pure state 1^) is entangled if it is not separable. 

A. Von Neumann entropy 

Pure bipartite entanglement is measured in terms of 
the reduced Von Neumann entropy S. For a pure state 
the reduction of its density matrix p = \ip)(ip\ is obtained 
through partial trace of one the two partitions as pa = 
trs/9 or, equivalently, ps = tr^p. Then, S is defined as 

S = S A = S B = -tx B {pB logps) ■ (1) 

The Von Neumann entropy provides an unambiguous 
measure of entanglement for a bipartite system in an 
overall pure state. For a separable state S = while 
for a maximally entangled state S — log JV, where J\f = 
min( A/a, TV's), with Ma = dim (7^) andA/s = dim(7Y B ). 
In what follows we will take the logarithm in Eq. base 
2. Thus, for a many-qubit system, the maximum value 
that the Von Neumann entropy can take is equal to the 
number of qubits that have not been traced out to obtain 
the reduced density matrix. 

B. Concurrence and entanglement of formation 

For the case of mixed states, the Von Neumann en- 
tropy is no longer a good measure of entanglement. If 
we consider a bipartite system on an overall mixed state, 
then each subsystem can have non-zero entropy even if 
there is not any entanglement .38). In order to mea- 
sure the bipartite entanglement of a mixed state we 
shall consider the so-called entanglement of formation 
Ep Starting from a mixed state with density ma- 

trix p = "YlnPityi) {">pi\, Ep{p) is defined as the average 
entanglement of the pure states of a given decomposi- 
tion of the mixed state p, minimized over all its possible 
decompositions: 

E F (p) = min . (2) 

{pi,ipi} 

Here E(\tpi)) is the amount of entanglement of the pure 
state measured, as discussed in the previous subsec- 
tion, by the reduced Von Neumann entropy. 
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Eq. J3J is operationally difficult to handle, since it in- 
volves an extremal condition. However, for the case of 
two-qubit systems, Ep can be expressed in terms of a 
much more amenable quantity, the so-called concurrence 
C E3. We have 



E F {p)=h[\\l + y\ -€■(,,)■- 



(3) 



where h(x) is the so-called binary entropy function de- 
fined by 



h{x) 



-xlog 2 a;- (1 - x)log 2 (l - x) 



(4) 



The concurrence C{p) of the two-qubit state p is defined 
as 



C(p) = max{0, c A } , 



(5) 



where c\ = X\ — A 2 — A3 — A4, the {Ai} being the square 
roots of the eigenvalues of the matrix p(a y (g> a y )p*(a y ® 
a y ), in decreasing order, and a y a Pauli matrix. Note 
that in this definition the complex conjugation is taken 
in the computational basis {|00), |01), |10), |11)}, where 
|0) and |1) are the eigenstates of <r z . From Eq. we 
see that Ep depends monotonously on the concurrence, 
which takes values between for separable states and 1 
for maximally entangled states. Moreover, it is easy to 
see that c\ take values in [—1/2, 1]. Thus, it is clear from 
Eq. © that a state is separable if c\ < and entangled 
otherwise. 

Other measures for pairwise entanglement exist. 
Among them we mention the entanglement of distilla- 
tion > the negativity 0] and the relative entropy EE] • 
All these quantities are related in one way or another to 
the concurrence. Therefore, we have chosen to present 
our results in terms of the concurrence. Nevertheless, we 
mention that we have also measured the negativity and 
found that it essentially gives the same results as the 
concurrence. 

Finally, for the case of multipartite entanglement the 
problem is much more subtle. Different measures of mul- 
tipartite entanglement have beenproposed, giving differ- 
ent results, even qualitatively [44|. Given the state of 
affairs, we will limit ourselves to discuss the existence 
of multipartite entanglement in terms of the qualitative 
information that can be extracted by comparing the av- 
eraged Von Neumann entropy for subsystems of different 
sizes. 



III. TRANSITION TO QUANTUM CHAOS 

Random Matrix Theory (RMT) was introduced to de- 
scribe the spectral properties of complex heavy nuclei. 
The key idea behind RMT is to replace the full physical 
description of the Hamiltonian by a suitable statistical 
representative of its symmetry group EE]- The statisti- 
cal spectral fluctuations of almost any complex Hamilto- 
nian were found to be described by a few classes of ran- 
dom matrix ensembles. This approach turned out to be 



very successful. The RMT analysis has been applied to 
many fields of physics such as nuclei, atoms, molecules, 
quantum dots, quantum billiards and many-body sys- 
tems among others O EE EE E3, EE H^l In the early 
1980's it was conjectured that the quantum versions of 
integrable and chaotic classical systems were described 
by different classes of random ensembles [5(1 l5l| . Since 
then, RMT has been successfully applied to describe the 
emergence of quantum signatures of chaos. 

The global manifestation of the onset of chaos in quan- 
tum systems consists of a very complex structure of the 
quantum states as well as in spectral fluctuations that 
are statistically described by RMT [45j . Let us focus on 
many-particle systems with two-body interaction as this 
is the nature of the model systems that we study in this 
paper. For this kind of systems it has been found that, 
under very general conditions, if the interaction strength 
exceeds some critical value, fast transition to chaos oc- 
curs in the Hilbert space of many-particle states. 

To be more precise let us consider a generic perturbed 
quantum many-body system. The Hamiltonian can be 
split into two parts: 



H = H o + V , 



(6) 



where Hq corresponds to the unperturbed original Hamil- 
tonian and the perturbation V to an interacting term. 
The unperturbed Hamiltonian Hq is assumed integrable. 
In other words, when V = the existence of as many inte- 
grals of motion as degrees of freedom is assumed. We will 
take the unperturbed eigenstates of the non-interacting 
Hamiltonian Hq to span the many-particle Hilbert space. 
In this basis, when the interaction is turned on, the 
eigenstates start to mix. The mixing can be described, 
for small interaction strengths, by perturbation theory. 
However, perturbation theory breaks down and quan- 
tum chaos sets in when the typical interaction matrix 
element between directly coupled states becomes of the 
order of their energy separation |l6j | (we say that two 
many-body states \ipi) and \ip2) are directly connected if 
(ipi\V\4>2) 7^ 0). The transition to chaos reflects in the 
statistical properties of the spectrum. 



A. Nearest neighbor level spacing distribution 

The nearest neighbor level spacing distribution P(s) is 
the probability density to find two adjacent levels at a 
distance s. 

For an integrable system the distribution P(s) has typ- 
ically a Poisson distribution 



Pp(s) = exp (-s) 



(7) 



In contrast, in the quantum chaos regime, for Hamilto- 
nians obeying time-reversal invariance, the nearest neigh- 
bor spacing distribution corresponds to the Gaussian Or- 
thogonal Ensemble of random matrices (GOE). This dis- 
tribution is well-approximated by the Wigner surmise, 
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which reads 

, . ITS ( TVS 2 \ , „ 

Pwd(s) = — expf — —\ . (8) 

Note that the GOE distribution exhibits the so-called 
"level repulsion" , i.e., the probability to find close energy 
levels is very small. This is in contrast to what is observed 
for integrable systems which exhibit level clustering. An 
example of Poisson and Wigner distributions is provided 
in Fig.El 

In Eqs. (JZJ) and (JHJ) the level spacing s is given in units 
of the mean level spacing A that we have set to 1. 

The transition to quantum chaos may be detected by 
the change of the nearest neighbor spacing distribution 
P(s) from Poisson to GOE. In order to obtain a more 
quantitative indication of this transition, it is useful to 
compute the parameter 

7 f° [P P (s)-P WD (s)]ds ' U 

where so ~ 0.4729 corresponds to the lowest s-value at 
which the Poisson [Eq. ((JJ] and Wigner [Eq. iJHJ] curves 
cross. This parameter takes values 1 and for the Poisson 
and Wigner distributions, respectively. 

The distribution P(s) describes the behavior of the 
fluctuations at energy scales of the order of A. Therefore, 
P(s) is a short range correlation. The effects of the onset 
of quantum chaos are also seen in higher moments of 
the distribution of energy levels. However, we will limit 
ourselves to the analysis of P(s) as a spectral signature 
of the transition to quantum chaos. 



B. Participation number 

The effects of the onset of chaos can also be observed 
in the eigenfunctions. The transition reflects in the de- 
gree of mixing of the eigenfunctions of the system. How- 
ever, the eigenfunction-based measures are more subtle. 
This is because the mixing of eigenfunctions is a basis- 
dependent quantity. Clearly, if the eigenfunctions are 
expanded in their own basis, they are not mixed at all, 
independently of the fact that the distribution P(s) is 
Poisson or GOE. Nevertheless, for Hamiltonians of the 
type © the increase of typical eigenfunctions' mixing 
with perturbation is naturally obtained in the basis {\i}} 
of the unperturbed Hamiltonian Hq (H^ . 

In the basis {\i)} the mixing (or, equivalently, the de- 
localization) of a given eigenstate 1-0) is customarily mea- 
sured in terms of the number of components (i\ip) which 
are significantly different from zero: A useful quantity to 
measure the degree of derealization of a given eigenstate 
is the so-called Participation Number (PN), defined as 



If the state is maximally localized, all its components 
are zero but one, with value Kil^)! 2 = 1 due to normal- 
ization. Therefore, in the localized regime £ « 1, while £ 
increases with increasing mixing. In the thermodynamic 
limit, £ is unbounded. However, in the case of finite size 
systems £ < N , where N is the dimension of the Hilbcrt 
space. For the GOE statistics, the PN is upper bounded 
by the value £ = N/3, due to the statistical properties of 
the chaotic states. 

Even though, in general, the transition from localized 
to delocalized eigenfunctions occurs in parallel with the 
transition from Poissonian to Wigner spectra, this is not 
always the case. As we will study in the following sec- 
tions, for particular model systems the spectra can be 
uncorrelated even in situations in which the eigenfunc- 
tions are delocalized. It is common to term this situation 
as weak chaos. 



C. Chaos and entanglement 

It is worthwhile to discuss what are the expectations 
for the entanglement content in the nearly integrable and 
fully chaotic situations. Let us consider a many-particle 
system with a Hamiltonian as in Eq. JBJ) and let Af denote 
the dimension of its Hilbert space. In a given basis {|n}} 
the density matrix for an eigenstate of the system, \tp) = 
J2n c n\ n ), writes as follows: 

Pnm = (n\lp)(ip\m) = c n c* m , (11) 

where c n = (n\ip) are the components of the eigenstate 
in the {\n}} basis. As discussed in the previous subsec- 
tion, the nearly integrable case corresponds to a situa- 
tion of weak interaction. This implies that in the ba- 
sis of the unperturbed Hamiltonian Hq the eigenfunc- 
tions are localized, i.e., c n rs S n _ n * for some n* € [1,A/]. 
Therefore, the density matrix will have all entries nearly 
equal to zero, except for the diagonal matrix clement 
Pn*n* = \c n * | 2 ~ 1- At the other extreme of strong inter- 
action, in which quantum chaos has set in, the eigenfunc- 
tions are fully extended. In this situation the eigenstates 
can be considered as random states with uniformly dis- 
tributed components with amplitudes c„ 1/yAf and 
random phases. In this case, the density matrix can be 
written as 

p«diag(l/AA,l/AA,...,l/AO + ft , (12) 

where is a Af x Af zero diagonal matrix with random 
complex matrix elements of amplitude « 1/Af. 

Suppose now that we partition the Hilbert space of 
the system into two parts with dimensions AT a and Afs, 
where NaNe = Af. The reduced density matrix pa is 
defined as follows: 

p A = Tr B/9 = ^2c nAnB c* n , AnB \n A ){n' A \ , (13) 

n B 
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where | n) = \n A nB)- Therefore, in the integrable case the 
reduced density matrix pa is a Ma x Ma matrix given by 

pA ~ diag (0, 0, . . . , 0, 1, 0, . . . , 0) . (14) 

In contrast, in the chaotic case 

p A w diag (1/A/" A , 1/A/U, . . . , I/AA4) + fi A , (15) 

where il^ is a zero diagonal matrix with matrix elements 
of O {\JNb /M) (sum of Mb terms of order 1/M with 
random phases). 

With Eqs. I|14l) and (|15l) in hand, it is easy to calculate 
the values for the reduced Von Neumann entropy and the 
concurrence. We have 



S A 



for integrable regime 
log(A/"^) for chaotic regime 



(16) 



A better estimate of Sa is obtained by considering the 
ensemble of random states according to the Haar mea- 
sure: S A « log(M A ) ~ M A /(2M B log e 2) 0. 

For the case of concurrence the dimension Ma — 4. If 
in Eqs. I|14|) - (|15fl we neglect the matrix £Ia we obtain 



for integrable regime 
-1/2 for chaotic regime 



(17) 



Thus, in both integrable and chaotic extremes, we obtain 
from Eq. (JSJ) that the concurrence is zero. 



IV. CHAOS AND ENTANGLEMENT IN SPIN 
CHAINS 

In this Section, we discuss bipartite and pairwise en- 
tanglement measures in two quantum lattice spin models 
in which a transition to chaos has been previously found 
and characterized. Both models have been proposed as 
suitable model for quantum computers. For the sake of 
completeness in this section we review the known prop- 
erties of these models, in particular the onset of quantum 
chaos. In parallel, we present new results concerning the 
behavior of quantum entanglement. In section TlV Al we 
focus, for a two-dimensional spin lattice, on the behav- 
ior of the concurrence as quantum chaos sets in. These 
results qualitatively agree with those presented in sec- 
tion llVBl for a family of one-dimensional spin models, for 
which we present a much deeper analysis of the behav- 
ior of the concurrence and of the Von Neumann entropy 
across the transition to quantum chaos. 



A. Two-dimensional spin lattice 

We consider consists of L spin-1/2 particles (qubits) 
placed on a two-dimensional square lattice in the pres- 
ence of an external static magnetic field directed along z. 



1 2 3 4 5 6 7 

JL/5 



FIG. 1: Level statistics parameter 7 as a function of the cou- 
pling parameter J (in units of S/L), for the model of Eq. I|18[l 
for a 3 x 3 lattice, calculated from the energy levels in the 
spectral band centered at — Ao averaged over 2000 random 
realizations of Si and Jij. The other parameters are Ao = 1 
and S = 0.09. 



Nearest-neighbor spins interact via Ising coupling with 
random strength. The Hamiltonian of the system is 



i<j 



o, cr- 



(18) 



The operators Oi are the standard Pauli matrices acting 
on the z-th qubit. The second sum in the Hamiltonian 
runs over nearest-neighbor spins and periodic boundary 
conditions are considered, r, corresponds to the energy 
separation between the states of the qubit i. Jij is the 
interaction strength between the qubits i and j. The 
parameters Tj and Jjj are randomly and uniformly dis- 
tributed in the intervals [A — 5/2, A + 5/2] and [— J, J] , 
respectively. This Hamiltonian was proposed as a model 
of isolated quantum computer with hardware imperfec- 
tions [To) . 

Here we focus on the case 5, J <C Ao, which corre- 
sponds to the situation where fluctuations induced by 
lattice imperfections are relatively weak. In this case, 
the unperturbed energy spectrum (Jij — 0) of Hamil- 
tonian l|18|) is composed by L + 1 well separated bands, 
with inter-band spacing 2Ao. Each band corresponds 
to states with a given number of spins "up" and spins 
"down". The highest density of states is obtained for 
the central energy band and therefore we expect that 
quantum chaos shows up first there. When interaction 
is turned on, a transition to chaos takes place. A value 
for the chaos border J c for this transition was given in 
[To|: J c cx S/L. This border was corroborated in [T^| . 
where the emergence of Fermi-Dirac thermalization in 
the chaotic regime was studied. A careful and detailed 
analysis of the transition to chaos for this model and its 
dependence on the size of the lattice has been taken in 
HO, [13 For the sake of comparison with the behavior 
of the entanglement measures we repeat some of these 



6 




FIG. 2: Participation number £ as a function of JL/8 for the 
model of Eq. I18H for a 3 x 3 lattice, with parameter values as 
in Fig. Q The PN was obtained from the central eigenfunc- 
tion of the band centered at — Ao averaged over 2000 random 
realizations. The dashed line corresponds to Nb/3. 



previous results. We consider a square 3x3 lattice. 
We note that for lattices composed of an odd number 
of qubits there is not a central energy band but instead, 
two central bands centered at ±Ao. In what follows, we 
will consider the states from the central band centered at 
-A Q . 

We have numerically diagonalized Hamiltonian (|18|l for 
different values of the interaction strength. To study the 
transition to chaos we have obtained the spectral statis- 
tics in terms of the nearest neighbor spacing distribu- 
tion -P(s) as well as the structure of the eigenfunctions in 
terms of the participation number £. We have restricted 
our calculations to the energies and eigenstates encoun- 
tered in the central negative band. For weak interactions 
the energy domain of this band is clearly visible and we 
keep the same domain even for stronger interaction where 
the band structure disappears. 

In Fig. ^ the parameter 7 as a function of the interac- 
tion strength is shown for Ao = 1 and S = 0.09. When 
the strength of the interaction J increases the P{s) dis- 
tribution smoothly changes from Poisson (7 = 1) towards 
GOE (7 = 0). Thus, increasing the interaction between 
the qubits a transition to quantum chaos occurs. In 
Fig. n we observe that the crossover from integrability 
to chaos takes place in the interval between JL/8 « 1 
and JL/8 ~ 5. 

At the same time, the eigenfunctions start to mix. For 
weak interactions the eigenfunctions are strongly local- 
ized: The number of components in the basis of the 
unperturbed Hamiltonian is of the order of one. At 
strong interactions {JL/8 3> 1) the eigenfunctions are 
extended, having a large number of non negligible com- 
ponents. This mixing of the eigenfunctions is shown in 
Fig. [21 in terms of the PN for the same parameter values 
as in Fig. ^ I n Fig- [21 we see that the PN smoothly 
changes from 1 (localized regime) to its upper bound 



FIG. 3: Concurrence C„ as a function of the coupling pa- 
rameter JL/8 for the model of Eq. I|18p for a 3 x 3 lattice, 
with parameter values as in Fig. C n was obtained from 
the central eigenfunction of the central negative band of the 
spectrum and averaged over all possible pairs of: a) nearest 
neighbors qubits (n = 1) and b) next to nearest neighbors 
qubits (n — 2). 



value of Nb/3 (chaotic regime), where Nf, corresponds to 
the number of eigenfunctions with energies in the central 
negative band. As we have discussed, the factor of 1/3 
arises due to the symmetries of the chaotic Hamiltonian 
that are described by the GOE. 

We now turn our attention to the entanglement mea- 
sures. We have calculated the concurrence between near- 
est (Ci) and next to nearest (C2) neighbor qubits. For 
this purpose we have drawn 2000 random realizations 
of Si and Jij and diagonalized Hamiltonian (|18|) . Us- 
ing for each realization only the central eigenfunction of 
the central negative band we have calculated the mean 
concurrence averaged over all possible nearest neighbor 
pairs of qubits. In Fig. [3] we show C\ (squares) and C2 
(circles), averaged over all the random realizations as a 
function of the strength of the interaction. In the ba- 
sis of the unperturbed Hamiltonian the concurrence is 
strictly zero. For very weak interactions {JL/8 <C 1), the 
concurrence remains small. At the other extreme, when 
the interaction is very strong {JL/8 1) and quan- 
tum chaos has set in, the concurrence is also small, as 
expected. Quite interestingly, the maximum of the C\ 
concurrence is for JL/8 ~ 1, that is, in the region in 
which the crossover from integrability to quantum chaos 
takes place. In Fig. [3] we can also compare the behav- 
ior of C\ with that of C2. The concurrence of next to 
nearest neighbor qubits is noticeable smaller than that 
of nearest neighbor qubits. This is not surprising as the 
Ising interaction in Eq. I|18|l couples only nearest neigh- 
bor qubits. Therefore, one should expect that quantum 
correlations between qubits are stronger for qubits that 
are close than for those farther away. However, we find 
that C2 is not negligible everywhere over the domain of 
J investigated, except for the integrable and chaotic ex- 
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FIG. 4: Color density map for the concurrence G\ calcu- 
lated for each individual eigenfunction (index in the vertical 
axis), averaged over 200 random realizations. Ciis plotted as 
a function of the coupling parameter JL/S for the model of 
Eq. I18H for a 3 x 3 lattice, with parameter values as in Fig.0 
In the inset, we show Ci averaged over the central 1/3 of the 
eigenstates of each individual band. The different curves cor- 
respond to the bands for states with: 4 (solid), 3 (dotted), 2 
(dashed), 1 (long-dashed) and (dot-dashed) spins "up". 



tremes, where also C\ goes to zero. Similarly to Ci, the 
concurrence Ci has its maximum value for JL/S w 1. It 
is interesting to point out that, similarly to what happens 
in the QPT in the Ising model 0, , the concurrences 
C\ and C2 exhibit their maximum values close to the 
value at which quantum chaos sets in. However, besides 
this similarity, there are other aspects in which the en- 
tanglement at the onset of chaos behaves in a different 
manner than for the case of a QPT. In section llVBI the 
qualitative differences of the behavior of entanglement in 
a QPT an for the onset of quantum chaos that we observe 
will be discussed. 

As we have discussed in section ITTTI the onset of quan- 
tum chaos occurs when the typical interaction matrix 
elements between directly coupled states becomes of the 
order of their mean level spacing. Thus, the onset of 
quantum chaos is expected to be observed first at the 
spectral energies at which the density of directly coupled 
states is larger. For the models discussed in this and 
next sections this happens at the center of the spectrum, 
i.e., in the central bands. Consequently, we focus our 
discussion on the eigenstates corresponding to the cen- 
tral energy bands. It is worthwhile mentioning that our 
choice is different from the studies of QPT for which the 
transition is a property of the ground state. For the case 
of a transition to quantum chaos we show in Fig. 21 the 
value of C\ as a function of the coupling parameter J, 
calculated for each eigenstates of a 3 x 3 lattice. The 
band structure can be clearly seen in the behavior of the 
concurrence. Inside each band, the concurrence grows 
from zero at J ~ to a maximum close to JL/5 = 1, 



after which it decays again to zero as J is further in- 
creased. On the other hand, for the states at the border 
of the bands, including the ground and the most energetic 
states, C\ behaves quite differently. This is due to the 
fact that these states do not mix significantly with the 
other states, as confirmed by our numerical data com- 
puting the participation number (data not shown). In 
the inset of Fig. 21 we show C\ , averaged over the central 
eigenstates of each band. We observe that, for all bands, 
C\ reaches its maximum value for J close to S/L. The 
fact that for the bands farther away from the center of 
the spectrum C\ decays slower to zero as J is increased 
is a finite size effect that should disappear at the thermo- 
dynamic limit. However, for the ground state C\ grows 
linearly. While we do not discard that at the thermo- 
dynamic limit even the concurrence calculated from the 
ground state will behave in a similar fashion as for any 
other typical state, Fig. 21 clearly indicates that for a fi- 
nite system this is not the case. 

B. One-dimensional spin chain 

In this section we discuss the behavior of bipartite 
and pairwise entanglement in a family of one-dimensional 
spin 1/2 chains. Due to its lower dimensionality these 
models will allow us to characterize the behavior of the 
entanglement across the transition to chaos in a deeper 
fashion than for the previous model. We shall find the 
same behavior for the concurrence than before. Nev- 
ertheless, with these models we are able to analyze its 
dependence on: the distance in the lattice between the 
partners, the range of the interaction and the size of the 
chain. Moreover, for one member of this family of models 
the chaos border does not coincide with the derealiza- 
tion border. This will give us the possibility to compare 
the behavior of the concurrence in a regime of weak and 
of hard chaos. 



1. Definition of the model. 

We consider a system consisting on a linear chain of 
L 1/2 interacting spins, subjected to a static transverse 
magnetic field (along z) and to a circularly polarized 
magnetic field rotating in the (x, y) plane with frequency 
v, B{t) = {B- L cos(vt + (p),-B ± sm(vt + tp),B z ) 0. In 
the coordinate system, which rotates around the z axis 
with frequency v, the Hamiltonian of this system can be 
written as 

ff= 4E {5 k a z k + r>(cos^ a* - shop a v k )} 

t-\ (19) 

k=l 

S k = uj k — v where u k is the frequency of the precession 
of the A;-th spin in the B z field, stands for the Rabi 
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frequency corresponding to the rotating field and J k ,k+i 
denotes the strength of the Ising interaction between the 
spins k and k + 1. The operators a k are the standard 
Pauli operators acting on the fc-th spin. In the following, 
we will take for simplicity <p — 7r/2 and consider that the 
static field B z has a constant gradient a along the chain 
such that S k — ak. Thus, the Hamiltonian takes the form 



k=l 



L-l 



x ^2 Jk.k+\o z k o 



fc+1 



(20) 



k=l 



We assume that for all k the inequality SI ^> Sf. holds. 
Open boundary conditions are taken. In [l3[ . this model 
was proposed as a possible candidate for experimental re- 
alization of quantum computation. The gradient of mag- 
netic field provides a labeling of qubits in terms of their 
Larmor frequencies. Thus, it allows for a way to address 
each qubit separately. 

It is worthwhile mentioning that, besides the differ- 
ent dimensionality, there is a more striking difference be- 
tween this and the previous model: The existence of a 
constant gradient in the magnetic field gives rise to a L- 
independent threshold for the onset of (weak) chaos. In 
|13| the transition to quantum chaos and its implications 
to quantum computation were explored. Here we want 
to discuss the behavior of entanglement in this model. 

In order to apply the approach discussed in section ITTT1 
it is convenient to represent Hamiltonian l|2UI) in the ba- 
sis in which it is diagonal for non-interacting spins. In 
this so-called effective field representation, the one-body 
unperturbed Hamiltonian Hq takes the form 



(21) 



fc=i 



and the interaction Hamiltonian V can be written as 
V = Vdi aK + Hand + K>ff, where 



L-l 



Vdiag = 

Vband = 

Voff = 

with 



1 x ^ 

9 Jk,k+l 



fc=l 
L-l 



b k bk+lCTk<Tk + l, 



a k a k+1 (j v k (j v k+l , 



(22) 



fe=i 

L-l 



r Jk ' k + 1 ( a k b k+i<J V k <T z k+1 + a k+1 b k a z k al +1 ) 



fc=i 



n 



ak 



(23) 



As before, the quantities Jk,k+i stand for the Ising inter- 
actions between nearest neighbor spins. In what follows, 
we will consider the interactions to be completely ran- 
dom, i.e. Jfe^/t+i = J£, where £ is a random number 
uniformly distributed in the interval [—1,1]. This model 




FIG. 5: Nearest neighbor level spacing distribution P(s) for 
the AA model (l c — L — l)of Eq. 12511 for a chain of 12 qubits 
and: a) J/J c = 0.35 and 6) J/J c = 15, calculated from the 
energy levels in the central band of the spectrum and averaged 
over 10 random realizations. The solid lines correspond to the 
Poisson (a) and Wigner surmise (b) theoretical distributions. 



is known as the NN model from its nearest neighbor char- 
acter. 

As in the two-dimensional model (|18f) . for the unper- 
turbed (J = 0) case the spectrum possesses a band struc- 
ture. Each band is characterized by a constant number 
n of qubits in the state |0) and L — n qubits in the state 
|1). When L is even, a central band (around zero) exists. 
It consist of the many-qubit states with L/2 spins "up" 
and L/2 spins "down". The number of these states is 
given by 



N h = 



(L/2)! (L/2)! 



(24) 



As discussed at the end of previous section, we will only 
consider the energy levels and energy eigenstates corre- 
sponding to the central band of the spectrum. 

When J > the potential term V mixes the states in- 
side each band and among different bands: In the basis of 
Hq, Vdiag is diagonal. Instead, Vband couples states which 
are either in the same band or in next to nearest bands. 
Voff couples states which are in nearest neighbor bands. 
The mixing of energy bands triggers the transition to 
chaos. For a relatively weak interaction the eigenstates 
(in the basis of Hq) are localized, while for stronger inter- 
action the number of components significantly different 
from zero increases. The transition from strongly local- 
ized to extended states occurs very fast with the increase 
of the interaction and sets in when the strength of the 
typical interaction is of the order of the mean level spac- 
ing between directly coupled many-body states [14|. In 
PJ| the value for the derealization border was found to 
be J c sa 4a 2 /fl. 

However, as it was shown in ^3|, the NN model is pe- 
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FIG. 6: Participation number £ as a function of the coupling 
parameter J for the AA model of Eq. 125ft . for a chain of size 
L = 12. The PN is averaged over all eigenfunctions in the 
central band of the spectrum and over 10 random realizations. 
The dashed horizontal lines correspond to Nt/3 and to N/3 
where TV = 2 L is the dimension of the Hilbert space. 



culiar in the following sense: The delocalization border 
does not coincide with the chaos border. Increasing the 
strength of the interaction J, the system goes from a 
regular regime to a weak chaos regime where the eigen- 
functions are delocalized but the level statistics is not yet 
described by Random Matrix Theory. If the interaction 
is further increased the bands overlap and the system en- 
ters a regime of strong chaos. This peculiarity is removed 
if the range of the interaction is larger than nearest neigh- 
bor. 

Here we will consider a range of the interaction l c from 
1 (for the NN model) up to L — 1 . The interaction term V 
keeps the same structure as in Eq. (122ft but the different 
terms are now 

L-l j + la<L 

Vdiag = JjkbjbkO-jCrl, 

3 = 1 k=j+l 

L-l j + la<L 

Vband = ~2^I X! JjkO-jaka^o-l, (25) 

j=l k=]+l 

- L-l 3+lc<L 

Voff = 2 XI H J jk(ajb k o-Vo-l.+akb j o-?o-l) . 

3=1 k=j+l 

For l c = L — 1 this model is known as the AA (All to All) 
model as in this case all qubits are allowed to interact 
with each others. In contrast with the NN model (l c — 
1), if l c > 1 the chaos border occurs at the same value 
J cr ~ 4a 2 /fl as the delocalization border [lflj . 

2. The onset of quantum chaos. 

In Fig.EJthe nearest level spacing distribution is shown 
for the AA model for a chain of 12 qubits. The P(s) dis- 



FIG. 7: Concurrence as a function of the coupling parameter 
J for the model of Eq. (12">> with interaction range l c — 5 for a 
chain of 10 qubits. For each eigenfunction C„ was obtained as 
the average concurrence between all possible pairs of qubits 
separated by a distance n: n = 1 for nearest neighbors pairs 
(diamonds), n = 2 for next to nearest neighbors pairs (right- 
triangles), n — 3 (circles), n — 4 (squares) and n = 5 (up- 
triangles) . The plotted values correspond to C n averaged over 
all the eigenstates of the central band and over 30 random 
realizations. 



tribution was obtained from the energy levels contained 
in the central band and averaged over 10 different ran- 
dom realizations. In panel a, the case of weak interac- 
tion (J/J c — 0.35) is shown. It is in good agreement 
with the Poisson distribution (solid line), as expected in 
the integrable regime. In contrast, panel b shows the 
situation corresponding to strong coupling (J/J c = 15) 
in which the P(s) distribution follows the Wigner sur- 
mise expected for a chaotic system with GOE statistics. 
In panel 6, the level repulsion effect is evident. Thus, 
when the interaction strength J > J c , the spectral level 
statistics changes from Poisson to GOE showing that a 
transition to quantum chaos is happening. 

Simultaneously, a localization-delocalization transition 
occurs for the eigenstates in the central band. This tran- 
sition takes place for any value of l c . However, for the 
NN model, the PN does saturate at a value which is lower 
than Nf,/3 corresponding to the case of Gaussian fluctu- 
ations. In Fig. El the PN is shown for the AA model for a 
chain composed of 12 qubits. The PN (squares) is aver- 
aged over all the eigenstates in the central band and over 
30 different random realizations. For weak interactions, 
the eigenstates are effectively localized (£ » 1). The 
PN increases monotonously with the interaction until it 
reaches the value N b /3 = 308 (lower dashed line). A com- 
plete mixing of different bands occurs for much stronger 
interactions (J/J c ~ 1000). This is seen in Fig. El where 
£ increases again and reaches its upper bound value cor- 
responding N/3, namely to one third of the dimension of 
the whole Hilbert space, as expected from RMT. 

The model of Eq. (|25|l shows a clear transition to quan- 
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FIG. 8: Mean concurrence C3 between qubits at a distance 
n = 3, averaged over all the eigenfunctions in the central 
band of the spectrum and over 30 random realizations. The 
mean concurrence is plotted as a function of the coupling 
parameter J for the model of Eq. 125H in a chain of 10 qubits. 
We compare C3 for different range of the interaction: The NN 
model corresponding to l c — 1 (diamonds), l c — 2 (triangles), 
and l c — 3 (circles). 



FIG. 9: Mean concurrence averaged over all possible pairs of 
qubits for all the eigenstates in the central band of the spec- 
trum and over 30 random realizations for the AA model of 
Eq. 1251 (l c = L — 1) as a function of the coupling parameter 
J and for different sizes of the chain of qubits: L = 6 (dia- 
monds), L — 8 (triangles), L = 10 (circles), L = 12 (squares). 



turn chaos in which both energy levels and eigenstates 
change their character. Now we turn our attention to 
the behavior of quantum entanglement. 



3. The concurrence: Sharing of entanglement. 

We have calculated the mean concurrence averaged 
over all the eigenstates in the central band as a func- 
tion of the strength J of the interaction. In Fig. [7| the 
mean concurrence is shown for the model of Eq. I|25l) for 
a chain of L = 10 qubits, with an interaction that couples 
L/2 neighbor qubits, i.e., l c = 5. The diamond symbols 
correspond to the mean concurrence C\ between nearest 
neighbor qubits . As it was discussed in section Ull CI we 
observe that C\ is close to zero in both extremes of chaos 
and of integrability. Moreover, similarly to what we have 
observed for the two-dimensional model l|18|) , in between 
the integrable and the chaotic extremes the concurrence 
increases and its maximum value is close to the value for 
the chaos border. Despite the fact that we have observed 
this behavior of concurrence for just two different models 
we conjecture that it is generic for transitions to chaos. 
In Fig. [7] the mean concurrence averaged over all qubits 
at further distances, C2 (right-triangles), C3 (circles), C4 
(squares) and C5 (up-triangles) are also shown. The be- 
havior of C\ and C2 is similar to that observed for the 
two-dimensional model. We see that the concurrence C n 
decreases with the distance n, except for weak interac- 
tions for which C2 > C\. This latter is a peculiarity of 
this model. It is interesting to notice that the coupling 
strength J max (^) at which the concurrence C n takes its 
maximum value does not change significantly with n. 



Let us now consider the following question: for a given 
distance n, how does C n varies if the range of the inter- 
action increases? In Fig. |H1 the mean concurrence C3 is 
plotted for different ranges l c of the interaction: l c = 1 
(diamonds), l c — 2 (triangles) and l c = 3 (circles). To 
this purpose, we computed C3 for all the eigenfunctions 
in the central band of the spectrum of a chain of size 
L = 10 and averaged over 30 different random realiza- 
tions. From Fig. |8] we observe that as the range of the 
interaction increases the mean concurrence C3 decreases. 
The same conclusions were also obtained for the behav- 
ior C n with n / 3 (data not shown). This fact can be 
understood from the pairwise character of the concur- 
rence. Since the amount of entanglement between one 
definite qubit and the rest of the system is bounded, this 
finite amount of entanglement has to be shared between 
all possible partners. When the range of the interaction 
is enlarged, it becomes easier for each qubit to become 
entangled with more qubits in the chain. As a result the 
pairwise entanglement between one single qubit and the 
rest of the chain is shared among more partners. There- 
fore, the average entanglement shared between two qubits 
decreases. This argument is valid if a change in the range 
of the interaction does not significantly change the total 
amount of bipartite entanglement that is shared between 
a single qubit and the rest of the chain. 

Finally, we have studied the behavior of the concur- 
rence as a function of the size of the system. In Fig. [5] 
the mean concurrence C a obtained from all possible pairs 
of qubits is shown for the AA model for different sizes of 
the chain: L = 6 (diamonds), L = 8 (triangles), L = 10 
(circles), L = 12 (squares). In this case, instead of mea- 
suring the concurrence C n for some value of n, we have 
measured the concurrence C a as the mean concurrence 
between all possible pairs of qubits in the chain. This 
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FIG. 10: Normalized distribution P(c\) of the quantity c\ — 
Ai — A2 — A :i - A 4 of Eq. JSJ for the AA model Eq. JSEJ with 
L = 12 and coupling: o) J/J c = 0.01 and 6) J/J c = 100.0. In 
the inset the behavior of the first moment of the distribution 
< ca > is shown as a function of the coupling parameter J. 
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FIG. 11: Von Neumann entropy Si between one qubit and 
the rest of the system; Si is averaged over all qubits and over 
all eigenstates in the central band of the spectrum and over 
30 random realizations for the AA model of Eq. 12511 as a 
function of the coupling parameter J and for different sizes 
of the chain of qubits: L = 6 (diamonds), L = 8 (triangles), 
L = 10 (circles), L = 12 (squares). 



is because, for the AA model, the concept of distance 
turns out to be meaningless, since the strength of the 
interaction between two qubits does not depend on their 
distance. The behavior of C a as a function of the size of 
the system can be well understood in terms of the same 
argument used to explain Fig. [SJ When the system size 
is increased, the number of possible partners with which 
a given qubit can be entangled also increases. Therefore, 
the pairwise entanglement decreases, in agreement with 
the data of Fig. We note that, for small system sizes 
(L = 6,8), the concurrence C a does not go to zero at 
the chaotic side of the transition. This finite-size effect 
disappears already for L — 10, 12. Similar results were 
obtained for the other models with different range of the 
interaction. 

The results presented in this section show that the 
concurrence maximizes for values of J which are close 
to those at which the transition to chaos occurs. As 
discussed in the previous section, the same behavior for 
the maximum of concurrence has also been observed for 
quantum phase transitions occurring in integrable mod- 
els (see, e.g., 0, ^ for a study on the Ising chain). 
However, as it can be seen in Fig. [S] for the transition 
to chaos, the concurrence approaches zero when the size 
of the system increases, in contrast to what is observed 
in QPT [2]]], where C n — > for all n > 1 but remains 
finite for n = 1. In these studies a critical scaling for 
the derivative of the concurrence was obtained. On one 
hand, the fact that for the onset of quantum chaos the 
concurrence diminishes when the system approaches the 
thermodynamic limit makes a finite size scaling analysis 
rather difficult, as numerical errors become soon of the 
same order of the measure itself. On the other hand, it 
is not clear whether the transition to chaos in the models 



we have considered becomes sharp at the thermodynamic 
limit. This means that to prove the existence of a criti- 
cal point for the transition remains an open problem. In 
order words, we lack for a critical point at which to per- 
form a scaling analysis. We have nevertheless analyzed 
the dependence of dC / d J on the system size without hav- 
ing found any clear indication of a scaling behavior (data 
not shown). 

It is interesting to study the different character of 
mixed pairwise entanglement in the integrable and 
chaotic sides of the transition. In Sec. lIII CI we gave sim- 
ple arguments that explain the different structure of the 
reduced density matrix p for a mixed bipartite state in 
the regimes of integrability and chaos. In the integrable 
region, due to the localized nature of the eigenstates, p 
is essentially diagonal with only one matrix element sig- 
nificantly different from zero. On the other hand, in the 
chaotic region, p is almost diagonal with matrix elements 
of comparable magnitude along the diagonal. Both cases 
give a very small (or zero) concurrence. However, while 
in the integrable case this is due to the fact that the 
two-qubit subsystem under investigation is essentially in 
a separable pure state, in the chaotic case the pairwise 
entanglement is zero due to the random structure of the 
eigenfunctions of the whole L-spin system. As a conse- 
quence, the two-qubit reduced density matrix is essen- 
tially diagonal. Therefore, in the chaotic regime, the 
interaction with the rest of the system mimics a deco- 
herence process for the two-qubit subsystem. 

The different origin of the very small value of concur- 
rence is illustrated by the distribution of the c\s (we 
remind the reader that the concurrence is defined as the 
maximum between c\ and 0, see Eq. JSJ). In Fig. 1101 
the probability distribution P{c\) is shown in: a) the 
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FIG. 12: Von Neumann entropy S n between left and right 
blocks of sizes n and L ~n qubits, respectively, averaged over 
all eigenfunctions in the central band of the spectrum and 
over 10 random realizations. The data correspond to the AA 
model of Eq. 1251 for a chain of size L = 12. From bottom to 
top, the different curves are from n = 1 to n = 5. In the inset 
S n (J I Jc = 15) is plotted as a function of n. The dashed line 
corresponds to a linear fit with a slope of 0.93 ± 0.02. 



FIG. 13: Von Neumann entropy as a function of the 

strength of the interaction J, averaged over all eigenfunctions 
in the central band of the spectrum and over 10 random real- 
izations. The data correspond to the AA model of Eq. 12511 for 
a chain of sizes L — 6 (diamonds), L = 8 (triangles), L — 10 
(circles), L — 12 (squares). In the inset Sl/2{J/Jc = 15) is 
plotted as a function of L. The dashed line corresponds to a 
linear fit with a slope of 0.52 ± 0.02. 



integrable regime and b) the chaotic regime, for the AA 
model and a chain of size L — 12. Clearly, in both cases 
the probability to find c\ > is very small. Therefore, 
the concurrence is very small in both cases. However, 
the distributions P(c\) are quite different. We note that 
numerical results about the distribution P(c\) in a differ- 
ent model of quantum chaos where presented in |35j | . In 
the Inset of Fig.^J one can see that the first moment of 
P{c\) changes across the chaos border. It is an interest- 
ing open problem to obtain an analytical form of P{c\) 
for integrable and chaotic situations. The possibility to 
use this distribution to mark the transition to chaos also 
deserves more investigation. 



4- The Von Neumann entropy. 

We now turn our attention to the behavior of bipartite 
entanglement measured in terms of the Von Neumann 
entropy First we consider the mean Von Neumann en- 
tropy Si of each qubit with the rest of the qubits in the 
chain. For this purpose, we divide the system in two par- 
ties: one consists of just one qubit and the other contains 
the remaining L — 1 qubits. Then, following Eq. lJT]l. we 
compute Si from the 2x2 reduced density matrix of the 
one qubit subsystem. In Fig.^Jthe behavior of Si across 
the transition to chaos is shown for the AA model and for 
different sizes of the system, from L = 6 to L = 12. We 
find that the bipartite entanglement Si shows the same 
behavior independently of the size of the system. The 
state of the system changes from separable to maximally 
entangled as the transition to chaos occurs. In all cases, 
the entropy saturates to its maximum value Si = 1, up 



to corrections of order 1/2 L |5^|. As discussed in Sec.lTTl 
these results show that there exists a global entangle- 
ment of each single qubit with the rest of the system 
and that this entanglement increases with the interac- 
tion. The maximum value of bipartite entanglement is 
obtained when quantum chaos has set in. 

We have obtained similar results for the bipartite en- 
tanglement when the two blocks in which the system is 
partitioned have different lengths. As an example we 
show in Fig. ll2l the Von Neumann entropy S n between a 
block of size n and the rest of the system (of size L — n) 
from n = 1 (bottom) to n = 5 (top) for the AA model 
and a chain of size L — 12. Similarly to the Si case, S n 
increases when the transition to chaos occurs and satu- 
rates to S n ~ n for large interaction strength. Hence, the 
state of the system becomes maximally entangled when 
chaos sets in. This is a direct consequence of the ex- 
istence of multipartite entanglement. Moreover, in the 
inset of Fig. ^1 we have plotted the saturation value of 
S n for J/J c = 15 as a function of the size of subsystems 
n. This shows that in the chaotic regime the bipartite 
entanglement scales linearly with the size of the smallest 
of the two blocks in which the global system has been 
partitioned: S n oc n. 

It is interesting to study the Von Neumann entropy 
as a function of the system size L, when the two blocks 
in which the system is partitioned have a size oc L. We 
have computed the bipartite entanglement Sl/2 corre- 
sponding to the case in which the system is partitioned 
in two halves. A value of Sl/2 > for any L is indicative 
of the existence of multipartite entanglement. The ob- 
tained behavior of St,/2 as a function of the strength of 
the interaction is shown in Fig. 1131 for different system 
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sizes. The behavior of Sl/2 is similar to that shown by 
S±. It takes very small values in the integrable regime 
and then increases with the interaction up to a value for 
which it saturates. The saturation value is w L/2 (up 
to corrections of 0(1) |53j). In the inset of Fig. we 
plot the value of Sl/2 for J/Jc = 15 (i.e., in the chaotic 
regime in which the eigenstates in the central band are ef- 
fectively mixed) as a function of the size L of the system. 
It is interesting to note that the Von Neumann entropy 
does feel the mixing of different spectral bands occurring 
for a very strong interaction (J/J c ~ 1000). The inter- 
band mixing (compare with Fig.|S]for L = 10) produces a 
increase in the Von Neumann entropy which nevertheless 
is small compared to that observed for the onset of chaos. 
This is in contrast with the pairwise measures such as the 
concurrence, for which we did not observed any change. 

In addition, for any given value of J when chaos has 
set in, the bipartite entanglement scales linearly with the 
size of the system. It is interesting to comment this result 
from the viewpoint of computational complexity. It was 
shown in [5S| that large entanglement of the quantum 
computer hardware is a necessary condition for exponen- 
tial speedup (with respect to classical computation) in 
quantum computation operating on pure states. To be 
more precise, a necessary condition for an exponential 
speedup is that the amount of entanglement increases 
greater than logarithmically with the size L of the com- 
putation. This condition is fulfilled in the chaotic regime 
where Sl/2 oc (L/2). We remark that, differently from 
problems like exact cover this is not limited to the 
transition region but extends to the whole chaotic regime. 
We also note that the relation between entanglement and 
computational complexity in quantum algorithms simu- 
lating quantum chaos has been investigated in Ref. |57|. 



C. Weak and hard chaos 

In this section, we discuss the behavior of quantum 
entanglement in situations of weak and of hard chaos. As 
it was discussed, the NN model of Eq. (|22|l . while similar 
in character to the model of Eq. I)25|l with l c > 1 , shows a 
quite unexpected peculiarity: The chaos border does not 
coincide with the delocalization border. Thus, when the 
strength of the interaction is increased, the NN model 
experiences a transition from integrability to a situation 
of weak chaos in which the eigenfunctions are delocalized 
while the level statistics is yet of Poissonian nature. This 
results from the fact that the NN Hamiltonian can be 
approximately mapped into a model of L free fermions 
as discussed in However, this non generic situation is 
removed if longer ranges of the interaction are considered. 
This gives us the possibility to compare the behavior of 
entanglement in situations of weak and hard chaos. 

We have calculated the nearest neighbor concurrence 
Cj for the NN model and for a long range interaction 
model with l c = L/2. In Fig. ^] we present our results. 
In panels (a) and (c) the 7 parameter and the PN £ are 




FIG. 14: Comparison between the NN model Eq. fHk (cir- 
cles) and the model of Eq. 1251 with l c — L/2 (squares), for 
a chain of qubits of length L = 10. ^From top to bottom, the 
three panels show, as a function of the coupling parameter J: 
the behavior of the spectral statistics parameter 7 (top panel); 
the mean concurrence Ci between nearest neighbor pairs of 
qubits (middle panel); the participation number £ (bottom 
panel). All quantities were calculated in the central band of 
the spectrum and averaged over 30 random realizations. 



shown, respectively. We can clearly see the peculiar be- 
havior of the NN model. For both NN (open circles) 
and l c = L/2 (solid squares) models the PN behaves in 
a similar way. For the NN model the PN signals a clear 
transition from localized to delocalized eigenstates, even 
thought it does not saturates at the value corresponding 
to Gaussian fluctuations. In particular, the delocaliza- 
tion border of both models coincide. Nevertheless, the 
7 parameter shows a very different behavior. For in- 
teraction strengths 0.2 < J/J c < 1 the level statistics 
parameter 7 for both models takes the value 7 « 1. This 
corresponds to the integrable case in which the nearest 
neighbor spacing distribution is Poissonian. However, for 
the NN model, 7 remains Poissonian for larger values of 
the interaction up to J / J c 100. This is the situation 
that has been termed as weak chaos. In contrast, for the 
l c = L/2 model the level statistics changes from Pois- 
sonian (7 = 1), to GOE (7 = 0), around J/J c = 1. 
Thus, for the l c = L/2 model the chaos border coincides 
with the delocalization border as it is commonly found 
in many-particle systems with two-body interaction. 

In panel (b) the corresponding results for the concur- 
rence C\ are shown. We observe again the difference in 
C\ due to the range of the interaction, as discussed in the 
previous subsection. Despite this difference, C\ shows a 
similar behavior for both models: It is small at both 
sides of the transition and increases in between, having 
its maximum value close to J / J c = 1. This numerical 
results suggest that the behavior of the pairwise entan- 
glement is more sensitive to the mixing of the eigenstates 
than to the onset of quantum chaos. 
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V. FINAL REMARKS 

We have studied the bipartite and pairwise entangle- 
ment in one and two dimensional spin lattice models that 
experience a transition to quantum chaos. 

To study the presence of multipartite entanglement, 
we have analyzed the behavior of the averaged Von Neu- 
mann entropy for subsystems of different sizes. In par- 
ticular, we have shown that, for a partition of the system 
into two equal-size subsystems, this quantity grows lin- 
early with the system size in the chaotic regime. This 
shows that the classical simulation method discussed in 
|55| cannot be used to efficiently simulate the quantum 
chaos regime on a classical computer. 

For the case of pairwise entanglement, we have studied 
the dependence of the concurrence on the distance be- 
tween the partners, the range of the interaction and the 
size of the system. Our results suggest that for a typical 
many-qubit state, the entanglement is mainly multipar- 
tite rather than pairwise. 

We have also discussed the different character that 
the pairwise entanglement has at the integrable and the 
chaotic side of the transition in terms of a suitable distri- 



bution of the eigenvalues of the two-spin reduced density 
matrix. The use of the moments of this distribution to 
mark the transition to quantum chaos remains and inter- 
esting open question. 

Finally, we have discussed the similarities and differ- 
ences between the behavior of the concurrence at a quan- 
tum phase transition and at the onset of quantum chaos. 
Our results show that the maximal concurrence is ob- 
tained close to the derealization border for which mixing 
of the noninteracting eigenfunctions takes place and not 
necessarily related to the onset of quantum chaos. 
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